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ABSTRACT 

We identify two classes of transiting planet, based on their equilibrium tem- 
peratures and Safronov numbers. We examine various possible explanations for 
the dichotomy. It may reflect the influence of planet or planetesimal scattering 
in determining when planetary migration stops. Another possibility is that some 
planets lose more mass to evaporation than others. If this evaporation process 
preferentially removes Helium from the planet, the consequent reduction in the 
mean molecular weight may explain why some planets have anomalously large 
radii. 

Subject headings: planetary systems: formation; binaries: eclipsing; scattering 



1. Introduction 

Ever since the discovery of short orbital period, Jupiter mass extrasolar planets (Mayor 
& Queloz 1995; Butler et al. 1997), the so-called 'Hot Jupiters', their nature and origin 
have been of great interest. This interest has only intensified with the discovery of transiting 
planets (Charbonneau et al. 2000; Henry et al. 2000), allowing the measurement of planetary 
radii and more accurate planetary masses. The identification of physical trends in the data 
may help to understand the physics of these objects. Some such trends recently identified 
include trends between planet mass (Mazeh, Zucker Sz Pont 2005) or gravity (Southworth, 
Wheatley & Sams 2007) and orbital period for the known transiting planets. The existence of 
such correlations presumably indicates something about the physical nature of these planets. 

In § [2] we investigate these correlations further, in order to identify the true physical 
variables that underlie them. We introduce a new set of planet interior models in §|2]for the 
purposes of comparing with the data. The resulting comparison to the observed trends is 
performed in § |H 
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2. Planetary Trends 



In order to understand the origins of the correlation between planet gravity or mass and 
orbital period (Figured]), we wish to examine the correlation further using variables that are 
more closely tied to the likely physics of these objects. Figure [2] shows the planet gravity 
plotted against the equilibrium temperature 



for the 17 known transiting planets with mass < 5Mj. This leads to an interesting difference 
with respect to Figure [1], in that the distribution shows a more bimodal character, with a 
hint of a gap at fixed Tgg between low and high gravity planets. To understand this further, 
we show in Figure [3] the same data, but now with gravity replaced by the Safronov number 



where Vesc is the escape velocity from the surface of the planet and Vorh is the orbital velocity 
of the planet about its host star. The division into two separate groups is now quite striking. 
There are two transiting planets that are discrepant. The planet HD 147506b is considerably 
more massive than the others and lies off the top of the plot in Figure [31 The isolated point 
at the lower left is the Neptune-mass companion to GJ 436b, whose measured radius is much 
below (Gillon et al. 2007) that of a H/He mixture, and may legitimately considered a different 
kind of object from the others. The rest of the planets define two separate classes, differing 
by almost a factor of two in their characteristic values of G at fixed Teq. We designate the 
planets with 6 ~ 0.07 ± 0.01 as Class I and those with 6 ~ 0.04 ± 0.01 as Class II. 

Inspection of the properties of the two classes reveals that the Class II objects orbit the 
hotter host stars in general, but at greater orbital separations than the Class I stars with 
similar Teq. This can be seen in Figure HI which compares the two parameters that determine 
the amount of planetary illumination. Class II also contains most of the planets who appear 
to have radii that are too large (e.g. Bodenheimer et al. 2001; Guillot & Showman 2002) 
compared to basic H/He structure model^. Nevertheless, the differences in planetary radii 
are not enough to explain the differences in B, because the radius anomalies are only ~ 10%m 
whereas the characteristic value of G differs by almost a factor of two between the two classes. 



^TrES-3 has the largest radius of the Class I planets. However, the slower cooling associated with the 
larger mass/heat capacity means that the observed radius agrees with the model at ages ~ 10^ years but not 
much beyond that. Thus, whether this planet can be considered anomalous or not is somewhat ambiguous. 




(1) 




(2) 



- 3 - 



The principal differences between tlie two classes are based on mass. Tlie planets of Class II 
are, on average, less massive than the corresponding planets in Class I, and orbit stars that 
are generally more massive. One way of illustrating this is shown in Figure O in which we 
plot the planet mass versus T^q. We see two different linear trends for the two classes, with 
the Class II planets being systematically less massive at equivalent Teg. Thus, not only is 
there a mass difference between the two classes, but there is a trend of planet mass with 
irradiation within each class. 



2.1. Formalising the two classes 

Setting aside then the Neptune-mass planet around GJ436b and the massive, eccentric 
planet around HD 147506 as different beasts, we are left with 16 planets, which we divide 
henceforth into two broad classes, based on the split in our Teg-6 diagram. The two classes, 
with the relevant parameters are given in Table [H The quantity T^g contains only stellar 
and orbit parameters, so the difference between the two classes in B is clearly determined by 
planetary properties, encoded in the ratio Mp/Rp. This fact is the principal empirical result 
that we wish to try and understand. The fact that the planets which have been identified as 
being anomalously large (e.g. HD209458b; OGLE-TR-10) fall mostly into Class II (although 
others, such as HD149026b, also fall in this class) is not responsible for the difference in G, 
but may offer a secondary indication as to the physical origin of this bimodality. 

In subsequent sections we want to try and determine the physical origin of this difference 
between Class I and Class II. To do so, first we introduce the models that we will use to 
compare to the data. 

3. Planet Models 

Our planet evolutionary models are based on the Henyey code of Hansen (1996; 1999), 
originally used to describe white dwarf evolution. Applying these models to the planetary 
evolution calculation is straightforward once one substitutes the correct input physics. The 
most important addition is the incorporation of the proper atmospheric boundary conditions 
for irradiated giant planets. For this we adopt the results of calculations using the PHOENIX 
code (Hauschildt & Baron 1999). The atmosphere models are used to generate a grid of 
boundary conditions ranging from effective temperatures of 2000 K to 50 K and surface 
gravities from 300 cm.s^^ to 3 x lO^cm.s^^. The grid was recalculated for each system, 
using the stellar host and orbital separation parameters in Table [TJ The irradiation from 
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the host star was included self-consistently in the radiative transfer in the manner described 
in Barman, Hauschildt & Allard (2001, 2005). The incident stellar flux was assumed to be 
distributed uniformly about the dayside hemisphere of the planet. Unless otherwise stated, 
the planet abundances were scaled to match the published metallicities of the host stars and 
molecular abundances were calculated assuming efficient rainout of dust (i.e. cloud free) 
as described in Barman et al. (2005). In most cases, rainout of dust resulted in severely 
reduced TiO and VO atmospheric number densities. Other important physical inputs are 
the use of the Saumon, Chabrier & Van Horn (1995) equation of state for H/He mixtures 
(already incorporated during the white dwarf applications) and the calculation of a detailed 
Rosseland Mean opacity table using the same opacities as used in the PHOENIX atmosphere 
models. At high densities, in the metallic phase, the H & He electrical conductivities were 
also included. 

Figure [7] shows an example of code verification, in which we reproduce the modelled 
evolution (dotted line) of the extrasolar planet HD209458b by Baraffe et al. (2003), using 
exactly the same boundary conditions (upper solid line). We also include a second model 
(lower solid line) which uses slightly different boundary conditions for the same system, 
calculated in the so-called 'rainout' approximation, rather than the 'cond' approximation 
used by Baraffe et al. The residual differences in the cooling curves likely reflect small 
differences in the Rosseland mean opacity tables used in the two cases. The good agreement 
shows that the theoretical models are consistent, although the gross discrepancy with the 
observed radius indicates that whatever underlying physics was missing from the Baraffe & 
Burrows models is also missing from ours. 

Using these models, we can repeat the exercises carried out by various other groups in 
the literature (e.g. Baraffe et al. 2005; Burrows et al. 2007) comparing our evolutionary 
models to the various planetary cases, using the proper stellar illumination to calculate the 
boundary conditions, in the rainout approximation. Performing the same calculations with 
our code yields similar results, so we will not reproduce all the results here. However, we 
do show, in Figures [8] and [9] the comparison between our models for eight speciflc planetary 
systems, four from Class I and four from Class II. In all cases, we calculate the models 
assuming that the planet has a cosmic composition of 73% Hydrogen and 25% Helium, 
using boundary conditions calculated with the correct illumination (using the appropriate 
separation and host stellar type for each object). In each case we include two curves. The 
flrst curve (dotted) indicates the radius calculated in the standard manner from the Henyey 
code. The second curve (solid) is calculated to include the additional 'transit radius effect' 
(Baraffe et al. 2003; Burrows et al. 2003), which takes into account that the surface of 
optical depth unity for limb-grazing stellar photons (which is what the transit observations 
actually measure) lies higher in the atmosphere than the location of the surface in our 
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model, which is taken to be the location where the Rosseland mean optical depth = 10^. 
Thus, we add to the model radius the vertical extent of the atmosphere between these two 
locations, taken from the Phoenix atmosphere models. The striking thing to note is that 
the models all match the Class I observations in that they either intersect the observed 
radius or lie above them (if the planet has a rocky core, the radius will be smaller than a 
coreless model of the same mass). The same is not true for the Class II planets, where the 
very same models that fit the Class I planets often fail to match the observed radii. This 
is the well-known anomaly that was first identified with the detection of the first transiting 
planet, HD 209458b (Bodenheimer et al. 2001; Guillot & Showman 2002). The fact that 
this phenomenon appears to be correlated with the difference between the two classes hints 
at a common physical origin for the anomalous radii and the difference in Safronov numbers. 



4. Possible Explanations 

We have identified a bimodality in the value of the characteristic Safronov number G 
for the bulk of the known transiting planets. This suggests that the physics that underlies 
the split is associated somehow with the planets ability to gravitionally scatter, capture or 
retain material. In this section we wish to examine how this might come about within our 
current understanding of planet formation, migration and evolution. 



4.1. Could it all be a selection effect? 

The cautious reader might wonder whether the fact that we find anomalously large 
planets around the hotter (and thus more massive stars) might be the result of some kind 
of systematic error in the estimation of the stellar radius. After all, the planetary radii 
in transit are measured relative to the stellar radius. However, these results come from a 
variety of groups and so it would require the error to be the same across several independent 
determinations i.e. it would require some fundamental community-wide misunderstanding 
of the structure of these stars. Figure [6] shows the ratio of radii versus the ratio of masses, 
i.e. the quantities that are actually determined by transit and radial velocity measurements. 
We see that there is a clear distinction between the two classes in this diagram, suggesting 
that the split cannot be resolved by simply rescaling and i?* for individual objects. 
Furthermore, the mass trend goes in the opposite direction. If the stellar radii were being 
biased systematically high, and then so would the stellar mass. Since planetary masses are 
also measured relative to the stellar mass, one would also expect them to be biased high. 
This is the opposite of the observed trend. 
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4.2. Energy Redistribution 

One of the parameters we have used to classify the stars into two groups is the so-called 
'equilibrium temperature Teq. This is important because, for giant planets so close to their 
host stars, the planetary structure (most importantly, the global entropy) is regulated by the 
irradiation it receives. To calculate this number we must make an assumption about how the 
incoming energy is redistributed across the surface of the planet. The calculation presented in 
equation ([1]) assumes that there is little redistribution over the surface of the planet, i.e. the 
area that is re-emitting with average temperature Teq is 2nRp. If the energy is redistributed 
efficiently across the surface of the planet, that area would be AirRp and the numerical value 
of Tf,q would be reduced by a factor ~ 0.84. Rescaling both classes by the same amount would 
obviously have no effect, but if one class of planet had efficient redistribution and the other 
did not, it could move the groups close together. There might even be some precedent, in that 
recent measurements of phase curves for a few extrasolar planets suggest different degrees of 
redistribution on planetary surfaces. Knutson et al. (2007) find that redistribution is quite 
efficient for the transiting planet HD189733b, a member of Class I. On the other hand, the 
measurement of a phase curve for the non-transiting planet v Andromedae (Harrington et 
al. 2006) suggests that redistribution is weak in the atmosphere of this planet. Although 
this system is not transiting and so not formally within our classification, the host star in 
this system is of spectral type F8 and thus plausibly a member of Class ij^. Furthermore, the 
large amplitude of the secondary eclipse in HD149026b (Harrington et al. 2007), a member 
of Class II, suggests little redistribution. 

However, scaling the Tf.q of one class would only result in horizontal motion of one of 
the groups in Figure [3] and would not explain the difference in O values. Even scaling T^q 
for Class II down by ~ 0.84 does not bring the two groups into alignment in Figure O Thus, 
differences in how the energy is redistributed cannot completely explain the difference in the 
two classes. Finally, the observations suggest that strong redistribution occurs in objects 
like HD 189733b, a Class I system, so that it is Class I that should have T^q scaled down, 
which would only increase the separation between the two classes. 



^ Since the inclination angle is unknown, we can place a lower limit on O for v And b by taking the largest 
plausible radius, yielding 6 > 0.045 for Rp = l.bRj. An upper limit may be obtained by taking a small 
radius and using the inclination constraint from the observed lightcurve {i > 30°). This yields 8 < 0.134 for 
Rp = IRj. Thus, formally we cannot rule out membership in either Class. 
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4.3. Evaporation 

The detection of excess Ly a absorption during primary transit of HD209458b (Vidal- 
Madjar et al. 2003) indicates that Hot Jupiter planets indeed lose some level of mass. 
However, the inferred lower limit on the mass loss rate is ~ VS^^g.s~^ ^ which is not enough to 
change the mass significantly. Before we discuss the details of possible evaporation scenarios, 
we first wish to explore whether evaporation can simultaneously explain the smaller masses 
and the larger radii of the Class H planets. The planetary radius is not particularly sensitive 
to mass, but what sensitivity there is results in an inverse relationship i.e. lower mass planets 
have larger radii than more massive planets under similar conditions. In particular, larger 
planets cool less rapidly, so one might imagine that the radius, which is essentially a measure 
of the planet's global entropy (e.g. Burrows & Liebert 1993), would be larger if the planet 
spent a significant fraction of it's lifetime with a higher mass than the present day mass. 
Can the lower masses of the Class II planets also explain their larger radii? 

Figure [TUl shows the radius evolution (solid hne) of a planet that has it's mass reduced 
at a constant rate from l.lMj to 0.62Mj over the course of 3 Gyr, while being subjected 
to illumination appropriate to HD 209458b. This mass reduction is of the order required to 
move an object from Class I to Class II. The dotted line corresponds to a planet that starts 
off at 0.62Mj and the dashed line corresponds to a planet where the mass loss is maintained 
at the same rate until the planet is completely eroded. We see that the evaporation leads to 
a slightly larger radius at intermediate ages, but that it only increases the planet radius by 
a few percent. It seems as though simple evaporation alone cannot explain the anomalous 
radii as a holdover from a prior, more massive and hotter state. 

4.4. Tidal Heating 

Another scenario that has been suggested for stopping the inward migration of planets 
is that the planet eventually overflows its Roche lobe (Trilling et al. 1998), possibly aided 
by tidal heating (Gu, Lin & Bodenheimer 2003). Could some process associated with this 
phenomenon give rise to the difference between the two classes? Indeed, Ford & Rasio (2006) 
note that the distribution of planetary mass with orbital period appears to be confined to 
lie outside a locus defined by twice the Roche limit for each system. This would arise if the 
final planetary configuration resulted from circularisation of an initially eccentric orbit, whose 
periastron was located at the Roche limit. Figure [11] shows the distribution of transiting 
planet mass ratios against semi-major axis. We see that, although the Ford & Rasio locus 
approximately delimits the observed sample, the trend is better fit by lines of constant 
Safronov number, as would be expected based on Figure [31 Certainly, the Roche limit does 
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not seem to offer any way to explain the difference between the two classes of planet. 

Could the inflation by tidal dissipation have an additional influence? In Figure [T2] we 
show the tidal luminosity required to inflate each planet to fill it's Roche Lobe at it's current 
position. This amounts to the assumption that the planet evolves outwards due to mass loss 
through the Lagrangian LI point, stopping at its present location when the planetary radius 
detaches from the Roche lobe. We see that there is some correlation between Lude and 6, as 
to be expected since Class I contains more massive planets on average, which in turn require 
more heating to inflate to a given radius. However, there is significant overlap between the 
two classes in terms of Lude, suggesting that this too is not the defining difference in the two 
classes. 



4.5. Accretion of Planetesimals 

It is perhaps illuminating that the split into two classes occurs when plotted in terms 
of the quantity 0, the Safronov number. This quantity is directly linked to the efficiency 
with which a planet gravitationally scatters other bodies (e.g. Safronov 1972). The most 
commonly accepted mechanism for planetary migration is loss of angular momentum due to 
torques in a gas disk (Lin, Bodenheimer & Richardson 1996). In the absence of such a gas 
disk, planets may still move inwards by either scattering other planets (Ford & Rasio 1996) 
or planetesimals (Murray et al. 1998). In particular, the planetesimal migration mechanism 
provides a natural stopping mechanism in that the planetesimals eventually collide with the 
planet rather than get ejected; a mechanism regulated by the Safronov number. 

Following Tremaine (1993), we can estimate the stopping criterion as follows. Scattering 
of planetesimals by a planet is a diffusive process, in which numerical integrations yield an 
energy diffusion coefficient per orbit (e.g. Duncan et al. 1987) 

^ i^^. (3) 

Consequently, the random walk in energy requires that a planetesimal undergo, on average, 
~ (l/aDE)^ periastron passages before being ejected. The probability that the planetes- 
imal strikes the planet on any given periastron passage is p ~ {Rp/aY / sini, where i is the 
inclination angle of the planetesimal orbit. Setting Np ~ 1 yields a value for the critical 
Safronov number 

• -x -3/4 

sm I ^ ' 



e,„i - 0.26 (^-^J . (4) 
Thus, O ~ 0.07 requires a planetesimal disk with an opening angle ~ 30°. Smaller G require 
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even larger opening angles, with a maximum 9 ~ 0.046, close to the characteristic value of 
Class II. 

This scenario offers several ways to get two classes of planet orbit. Perhaps the simplest 
possibility is that one class migrates via gas torques and another by scattering. If all planets 
migrate by scattering, then the separation in 6 requires that the population of planetesimals 
scattered by Class II planets be dynamically 'hotter' than those for Class I planets. Alter- 
natively, a two planet system scattering planetesimals can result in an effective repulsion 
between the two planets (this is thought to explain the outward migration of Uranus and 
Neptune in our own solar system e.g. Fernandez & Ip 1984). Perhaps the Class II planets 
have a second body in the system, that results in them being pushed in further. Another 
possibility is that the difference results from the Class I planets accreting a larger fraction of 
their mass in heavy elements as the migration slowed due to their inability to eject planetesi- 
mals from deep within the potential well. Finally, we note that the above diffusion coefficient 
comes from calculations performed for Oort cloud applications, where planetesimal accretion 
is a much smaller effect than it is likely to be in this case (Hansen 2000). The likelihood of 
these various scenarios can be better estimated once this calculation has been revisited. 

If the Class I planets accrete a significant amount of rocky material, that may offer an 
explanation for why some radii, but not all, are anomalous. Most models for the anomalous 
radius of HD209458b and other extrasolar planets postulate an extra source of energy of 
some kind (e.g. Lin, Bodenheimer & Mardling 2001; Showman & Guillot 2002; Winn & 
Holman 2005) or a way to retard the cooling more than normal (Chabrier & Baraffe 2007). 
If we postulate that such an energy source is present in all planets, we can explain the 
now smaller radii of the Class I planets by virtue of the fact that they have significantly 
supersolar metallicities, which means that their mean molecular weights are larger and their 
radii correspondingly smallejfl. This is qualitatively similar to the proposal made by Guillot 
et al. (2006), in which they postulate that a fixed fraction (0.5%) of the irradiance energy 
is somehow 'recycled' as internal energy, albeit by an unknown mechanism. 

To examine this in more detail, we recalculate our models including an additional heating 
term. The extra energy is deposited throughout the planet, proportional to local mass 
fraction. To illustrate the resulting effect, we recalculate the evolution of HD209458b. In 
order to match the observed radius, we need to postulate an internal luminosity ~ 2.3 x 
lO^^ergs.s"^, which amounts to 1.7 x 10~^ of the absorbed stellar flux. We can do a similar 



•^Burrows et al. 2007 have recently suggested that increased metallicity may make planets larger, but 
they only increased the atmospheric opacity, without increasing the mean molecular weight in the envelope. 
Inclusion of this latter effect, necessary for a self-consistent model, acts to make the radius of the planet 
smaller, not larger, as one increases the metallicity. 
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calculation for each of the transiting planets. There are now two adjustable parameters 
(once one assumes the mass and irradiation based on the observations) - the amount of 
extra heating and the core mass. The first variable makes the planet more extended, while 
the second makes the planet more compact. Thus, one can trade off the two against one 
another in order to fit the observed radius. Figure [13] demonstrates this for four representative 
planets in our sample. Each curve represents the relationship between the extra heating e 
and the core mass (in earth masses), determined by matching the observed radius of the 
planet. Here e is defined as a fraction of the irradiance for each planet, in the spirit of 
Guillot et al. (2006). In some cases, e at finite core mass. This occurs when the 
observed planet radius is smaller than the expected radius for a coreless model with no extra 
heating. On the other hand, some models also require a finite e even at zero core mass. This 
occurs when the observed radius is anomalously large. All models in Figure [T3] are calculated 
in the asymptotic limit (the extra heating contribution sets a floor in the entropy and hence 
radius below which the planet does not contract). Thus, TrES-2 has a non-zero e at fixed 
mass even though we did not classify it as anomalously large based on it's current age limit 



Examination of Figure [12] allows us to assess the viability of this model for explaining 
the observed difference in Safronov number/Mass. If we follow the Guillot et al. notion of 
a fixed irradiance fraction e, we can read off the core mass required. For example, a value 
of e ~ 3 ± 2 X 10^^, implies a core mass for TrES-2 of 33 ± 13M^, which is only ~ 10% of 
the total planet mass. Thus, in order for TrES-2 to have accreted almost 50% of its mass 
in heavy elements (as required to move it from Group II to Group I), the internal energy 
source would have to be much stronger than is usually assumed. A similar argument can 
be made for other Class I sources, such as TR-113, which requires a core mass 98 ± ISM®, 
which is ~ 23% of the total. In summary, the scenario in which the Class I planets accrete 
a large amount of planetesimals places even greater demands on the unknown energy source 
than in the traditional case. 

Could the accreted material supply it's own heat, via radioactive decay? Using the 
estimated energy output and abundance of ^^K for the Earth's core (Gessmann & Wood 
2002 and references therein), one can estimate an energy input of 



which is orders of magnitude too small to affect the structure of these planets. 

Thus, while it is certainly possible to make a consistent model in which all planets have 
a mysterious energy source but the Class II objects have accreted a substantial amount of 
material to reduce their radii, the results are unsatisfactory on two counts. One is that it 



of > 1 Gyr. 
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still requires the existence of a mysterious energy source, and the other is that the amount 
of mass accretion required for the Class I is prohibitive, requiring in turn a much stronger 
energy source than is invoked for the Class II sources. 



The model of § 14.51 invokes enhanced metallicity to reduce a planetary radius. The 
radius is reduced whether the metals are concentrated in the core or spread throughout 
the envelope. This occurs because the increase in metallicity increases the mean molecular 
weight of the envelope (or core) and leads to a steeper radial pressure gradient and a more 
compact planet. It then follows that decreasing the mean molecular weight can increase the 
planetary radius. Is it possible to increase the radii of Class II objects by reducing their mean 
molecular weight? Since the default mixture is one composed of 73% Hydrogen and 25% 
Helium, is it possible to enhance the Hydrogen fraction of the envelope and thereby increase 
the radius? For objects in the partially degenerate regime, such as planets, the radius scales 
R ~ SO that it should be possible to extract a radius increase of approximately 20% 
by replacing the Helium in the model planet with Hydrogen. This is more than required to 
explain the radius anomalies of the Class II planets (which are ~ 10%). 

As an example. Figure [T3] shows the effect of removing the Helium on the evolution of 
HD209458b. We use a boundary condition for an atmosphere computed assuming a solar 
(2%) metallicity but depleted in Helium. We also reduce the Helium content of the envelope 
by the same fraction. We see that global He mass fractions Y ~ 0.14±0.02 are consistent with 
the observed radius. In order to fit the other anomalously large planets, HAT-P-1 requires 
a model with Yf < 0.14 and for WASP-1, Yf < 0.13. However, not all Class II planets 
require low Helium content. For instance, TR-56, despite it's large radius, is well fit by the 
cosmic composition model, because of it's larger transit radius increment (a consequence of 
the higher surface irradiation). 

Thus, we have the possibility of solving both the anomalous radius problem and explain- 
ing the difference between the two classes, if we invoke significant evaporation in which the 
material that is lost carries a supercosmic abundance of Helium. To quantify this, consider 
the following simple model. We start with a planet that has initial mass Mq, with the correct 
cosmic proportion of Helium, Yq = 0.25. If we remove a certain fraction of Hydrogen and 
Helium in proportion a = AMh/ ^Mne, we are left with a final mass Mf and final Helium 
mass fraction Yf. The initial and final masses are related by 



4.6. Evaporation and Helium-poor Planets 




(6) 
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We know Mj and determine Yf by fitting our evolutionary models. If we furthermore 
require that Mq ~ 1.8Mj, in order to match the difference in mass between Class I and 
Class II, then we can infer the require value of a. For HD209458b, with Yj ~ 0.14, this 
yields a ~ 1.58, instead of the value a = 3 which would preserve cosmic composition. This 
corresponds to a wind in which Helium is 39% by mass or 14% by number. 

Figure [15] shows the effect of Helium-rich evaporation in an evolutionary model. We 
start with a IMj planet irradiated in the manner of HD209458b, and remove mass at a fixed 
rate until the planet reaches 0.63Mj, after which mass loss is shut off. The mass removed is 
composed of 50% Helium by mass, so that the global Helium mass fraction of the planet is 
reduced, reaching Y = 0.094 by the of the mass loss. We also show the effect of evaporation 
assuming a solar composition wind. The planet actually expands as Y drops. This expansion 
allows us to match the observations of HD 209458b. 

The final point to note is that, if the observed planet population were composed of 
systems in varying stages of evaporation, the distribution in would exhibit a continuous 
distribution, rather than the observed bimodality. In order to match the observations, the 
postulated evaporation must only last for a finite time, shutting off when ~ 0.04. 

5. Discussion 

We have identified a curious bimodality in the distribution of transiting planet proper- 
ties, in that they (apart from two outliers which may be further distinguished on the basis 
of other characteristics) seem to possess one of two distinct values of the Safronov number. 
We have explored several reasons for the dichotomy, two of which seem the most promising. 

The first possibility is that the Safronov number, which essentially measures the effi- 
ciency with which a planet scatters other bodies, plays an important role in determining 
when a planet halts its migration. This would emerge quite naturally in the case where some 
or all of the planets migrate as a result of ejecting smaller bodies through repeated grav- 
itational encounters. While this scenario naturally explains the prominence of G, it offers 
no convincing explanation for a secondary observation, namely that those transiting planets 
with anomalously large radii seem to fall predominantly into Class II, the class with smaller 
values of 9. 

It is interesting that our split into Class I and Class II planets show some similarities 
with the assertion of Burkert & Ida (2007), that the distribution of planets in semi-major 
axis interior to 1 Au is different for planets around the highest mass (> I.2M0) stars in 
the sample and that the planets in those systems are also systematically of lower mass. 
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These two features are also defining characteristics of our Class II planets, suggesting that 
the distinction carries over to greater separations than is currently probed by the transiting 
planet sample. This might be interpreted as giving support to the suggestion that multiple 
migration mechanisms are at work. 

A second possibility is that the difference in is a consequence of the fact that some 
of the planets have lost a markedly larger fraction of their mass through some form of 
evaporation. This scenario requires that many of the Class II planets lose ~ 50-60% of their 
initial mass (but not much more). The fact that Class II planets are found preferentially 
among those with hotter host stars is consistent with this idea, in that planet evaporation 
models are quite sensitive to the amount of EUV radiation received (e.g. hammer et al. 
2007 and references therein). The exact amount of evaporation is still difficult to calculate 
because it depends sensitively on the amount of chromospheric emission from the host star, 
a quantity poorly known even for F and G stars. This scenario also offers a way of explaining 
the anomalous radii, if we allow for the possibility that the evaporation reduces the Helium 
abundance in the planet, i.e. the mass in the wind carries a larger than solar abundance of 
Helium. The consequent reduction in the mean molecular mass of the planet results in a 
larger radius for a fixed entropy. The anomalously large planets can all be fit with Helium 
mass abundances Y ~ 0.14. We note that this need not happen in all cases, since some of 
the Class II planets can be fit with normal abundance models. 

How might such an evaporation occur? Inspection of the exosphere models of Yelle 
(2004) indicates that the evaporating wind from an irradiated planet consists of ionized 
Hydrogen but neutral Helium (for the exosphere temperatures ~ lO^K lie between the ion- 
ization temperatures of the two atoms). Thus, if a planet possesses a magnetic field, the loss 
of hydrogen may be reduced (although not completely cut off - there is an atomic Hydrogen 
component to the wind loss as well) in such a way that the composition of the mass lost may 
be enriched in Helium. An alternative model is that Coronal Mass Ejections from the star 
may have a significant impact on the mass loss (Khodachenko et al. 2006), again depending 
on the strength of the planetary magnetic field. Whichever process is in operation must also 
switch off when it reaches a limiting state defined by our Class II, since the observed dis- 
tribution is not continuous. Thus, we expect a somewhat different evaporation history than 
that outlined in Baraffe et al. (2005), in which planets evaporate down to Neptune masses. 

Where do Jupiter and Saturn fit into this picture? Figure [16] shows the evolution of 
Teg and 9 in the event that Jupiter or Saturn migrated closer to the Sun. We see that 
they could easily match the position of either group, depending on their location. This 
also illustrates the peculiarity of the gap between the two groups, as it would be perfectly 
reasonable to place Jupiter within the gap, if it was located at distances ~ 0.027 Au. On 
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the other hand, it is intriguing that, if Jupiter and Saturn were placed independently at 
the same close distance from the Sun, the difference between their values of would be 
reminiscent of the gap between Class I and Class II. Furthermore, it has long been claimed 
that Saturn's atmosphere is depleted in Helium, although the exact value is still somewhat 
uncertain (Gautier et al. 2006). For many years this has been understood in terms of the 
separation of Helium in the metallic Hydrogen core of the planet (Stevenson & Salpeter 1977) 
in order to explain Saturn's anomalous luminosity (e.g. Pollack et al 1977). Could similar 
physics be responsible for the difference between Class I and Class II? Fortney & Hubbard 
(2003) have already investigated the role of Helium phase separation in irradiated extrasolar 
planets and have concluded that the interior temperatures for the Hot Jupiters are likely 
too high for such processes to be relevant. It therefore appears that the Helium paucity of 
Saturn's atmosphere and it's likely membership in Class II had it migrated inwards, may 
simply be a coincidence. 

The last couple of years have seen a variety of measurements concerning the properties 
of these transiting planets, many of which probe the state of the planet atmosphere directly. 
The extreme difficulty of these measurements means that comparison to atmosphere models 
are essential and thus the interpretations are necessarily linked with what one assumes for 
the model. We hope that some of the above considerations will encourage modellers to take a 
very broad view of their input assumptions and cover the widest possible range of parameter 
space, including both heavily metal enriched atmospheres and atmospheres with strongly 
reduced Helium abundance. 
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Table 1. The properties of tlie known transiting planets, based on the tabulation of 
http://obswww.unige.ch/~pont/TRANSITS.htin, as of June 1 2007. Ages are taken from 
the compilation of Melo et al. (2006). Where no data was available, we have assumed a 

lower limit of 1 Gyr for the age. 



Planet Name 


Mp 


Rp 




Teff 


a 


Age 


T 

eg 


e 




(Mj) 


(Rj) 


(Mo) 


(K) 


(AU) 


(Gyr) 


(K) 






Class I 


OGLE-TR-Ul 


0.52 


1.01 


0.81 


5044 


0.047 


> 1.1 


1027 


0.062 


OGLE-TR-113 


1.35 


1.09 


0.77 


4804 


0.023 


> 0.7 


1345 


0.076 


HD189733b 


1.15 


1.16 


0.82 


5050 


0.031 


> 0.6 


1201 


0.079 


TrES-1 


0.76 


1.08 


0.89 


5250 


0.039 


2.5±0.1 


1151 


0.065 


TrES-2 


1.20 


1.22 


0.98 


5850 


0.037 


> 1 


1474 


0.077 


XO-1 


0.90 


1.18 


1.0 


5750 


0.049 


> 1 


1210 


0.078 


WASP-2 


0.88 


1.04 


0.79 


5200 


0.031 


> 1 


1292 


0.069 


TrES-3 


1.92 


1.30 


0.90 


5720 


0.023 


> 1 


1645 


0.078 


Class II 


OGLE-TR-10 


0.61 


1.22 


1.10 


6075 


0.042 


> 1.1 


1535 


0.040 


OGLE-TR-56 


1.29 


1.30 


1.17 


6119 


0.023 


3± 1 


2262 


0.040 


OGLE-TR-132 


1.14 


1.18 


1.26 


6210 


0.030 


> 1 


2007 


0.048 


HD 149026b 


0.33 


0.73 


1.3 


6147 


0.042 


2 ±0.8 


1743 


0.031 


HD209458b 


0.66 


1.32 


1.10 


6117 


0.047 


4.5 


1445 


0.045 


HAT-P-1 


0.53 


1.36 


1.12 


5975 


0.055 


> 1 


1318 


0.040 


WASP-1 


0.87 


1.44 


1.15 


6110 


0.038 


> 1 


1819 


0.042 


XO-2 


0.57 


0.97 


0.98 


5340 


0.037 


> 1 


1316 


0.046 



Unclassified 



HD147506 8.17 1.18 1.35 6290 0.069 > 1 1556 0.737 

GJ436 0.07 0.35 0.44 3200 0.028 > 1 612 0.027 
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Fig. 1. — The upper panel shows the planet gravity plotted against orbital period, for the 
known transiting planets. The lower panel shows the planet mass versus orbital period. In 
both trend is visible, with shorter period planets having higher masses & gravities 

on average. 
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T 




Fig. 2. — The planet gravity plotted against equilibrium temperature, as defined in the 
text. Note that the distribution is broader than when the gravity is plotted against period 
and that there appears to be a hint of a gap between high and low gravity planets. 
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Fig. 3. — If we replace the gravity with the Safronov number O, we find that there are now 
two clear groups at fixed equilibrium temperature T^g, apart from two outliers discussed in 
the text. We label them as 'Class I' (open points) and 'Class IP (solid points). 
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Fig. 4. — One might term this the 'Illumination diagram', as the two quantities plotted 
(stellar effective temperature and semi-major axis scaled by the stellar radius) determine 
how much stellar ffux is received at the surface of the planet. Open and filled circles indicate 
Class 1 and Class 11 planets, as in Figure [31 To plot a planet on this diagram one does not 
require the planetary mass or radius and so we can also plot the positions of all the other 
radial velocity planets as well. These are shown by the star symbols. Also shown are dotted 
lines indicating lines of constant T^q. The split into Class 1 and Class 11 is not as obvious in 
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this diagram, although the general trend is for Class II to orbit hotter stars. 




Fig. 5. — Planet masses plotted against equilibrium temperatures. The two classes of planet 
both appear to obey an approximately linear relation with T^q, but with slopes different by 
almost a factor of 3. 
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Fig. 6. — The open circles are Class I and filled circles are Class II. Rescaling an individual 
star will not move a point in this diagram, suggesting that the observed dichotomy is real. 
The absence of objects near (9 x 10^^,0.113) is curious, and unexplained. If this were 
indicative of some systematic error in the measure of Rp/ it would still not explain the 
dichotomy, since increasing the radius of the Class II planets would just make © smaller and 
similarly, the © of Class I would just get larger. 
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Fig. 7. — The dotted curve is the model of Baraffe et aL (2003), while the dashed curve is 
from Burrows et al. (2007). The upper solid curve is our model calculated using the same 
boundary conditions as Baraffe et al. (termed the 'cond' approximation). The lower solid 
curve is our model for the same planet, but with the boundary conditions calculated in the 
"rainout" (cloud-free) approximation, which we will use throughout the rest of the paper. 
Planet mass is assumed to be 0.69Mj in our models, as in the Baraffe model (the Burrows 
model is for a 0.64Mj planet). The observed value of the planetary radius is shown in the 
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upper right. Note that we have not included the transit radius correction in this plot, since 
the principal goal is the comparison of Henyey model results. 
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Fig. 8. — The four panels show evolutionary curves for four of the Class I planets, using 
masses and boundary conditions appropriate to each system. The dotted line indicates the 
traditional radius and the solid curve indicates the expected optical transit radius. The 
measured planetary radii are also shown, and are consistent with the models in all cases, as 
long as moderate heavy element cores are allowed. 
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Fig. 9. — The four panels show evolutionary curves for four of the Class II planets, using 
masses and boundary conditions appropriate to each situation. The dotted line indicates 
the traditional radius and the solid curve indicates the expected optical transit radius. The 
measured planetary radii are also shown. Unlike the Class I planets, several of these planets 
seem to be too large compared to the models, unless the planets are implausibly young. Note 
also that the transit correction can vary from system to system, as it depends somewhat on 
the level of illumination. 
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Fig. 10. — The solid curve shows the radius evolution of a planet evaporating from l.lMj 
to 0.62Mj over 3 Gyr. The dotted line indicates the cooling curve for a 0.62Mj planet. 
The boundary conditions used are appropriate for the HD209458b case and the filled circle 
indicates the observed value for that planet. 
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Fig. 11. — Filled and open circles once again indicate Class I and Class II planets. For the 
eccentric transiting planet HD 147506, we include an error bar to indicate the radial excursion 
between semi-major axis and periastron. The dashed lines indicate the Roche limit and twice 
the Roche limit. The dotted lines indicate lines of constant Safronov number (0 = 0.04 for 
the lower curve, © = 0.07 for the upper curve) assuming a planetary radius of 1.2 Rj. 
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Fig. 12. — The open and solid points again indicate Class I and Class II planets respec- 
tively. Ltide is calculated as the internal luminosity required to inflate a H/He planet of the 
appropriate mass to fill the Roche lobe of each planet at it's current location. (Note we have 
not plotted a point for HD149026b, since that is clearly significantly different from a H/He 
planet.) 
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Fig. 13. — The four curves show the relation between insolation recycling efficiency e and 
core mass Mcore^ as determined for four different planetary systems. In each case, the values 
are constrained so that the combination results in a model radius that matches the observed 
radius for that system. 
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Fig. 14. — The dotted line is the evolution of a 0.65Mj planet with a solar composition 
atmosphere and envelope, irradiated in a manner appropriate to HD209458b. The solid 
lines indicate models in which a fraction of the Helium has been removed from the envelope, 
leaving the planet with a Helium mass fraction Y as indicated. 
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Fig. 15. — The solid line shows the evolution of a planet whose mass is reduced from 1 Mj 
to 0.63 Mj over the course of 3 Gyr. The mass removed is 50% Helium by mass, so that 
the mean molecular weight of the planet drops while mass is being lost. The dotted line 
shows the evolution if the same amount of mass is lost over 1 Gyr. Finally, the dashed line 
shows the evolution if the mass loss occurs over 3 Gyr, but the mass removed has the normal 
cosmic composition of Helium. 
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Fig. 16. — The dotted and dashed curves show how Jupiter and Saturn would move in the 
diagram as one moved them closer to the sun. We have not adjusted the radii to 
account for the insolation as it is a small effect on 0. 



